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The evolution of an inhomogeneous universe composed entirely of matter is followed from an early, nearly 
uniform state until the time when the inhomogeneities have begun to grow large. The particular distribution of 
matter studied in this article is chosen to have a periodic variation in only one of the directions, which is simple 
enough to allow the behavior of the metric to be solved analytically, well beyond a linear approximation based 
on the initial smallness of the fluctuation. 

This example provides an illustration of a universe where the inhomogeneities can affect its average expansion 
rate; and its simplicity allows a condition to be derived that tells when their presence should begin to become 
important. Since the averages of the non-uniform parts of the metric and the matter density grow faster than their 
uniform parts, the average expansion rate accelerates with the advent of the era governed by the inhomogeneities. 

PACS numbers: 04.25.-g,98.65.Dx 



I. INTRODUCTION 

One very apparent feature of the universe is its nonunifor- 
mity, which is evident even over quite large distances. At the 
scale of a few hundred megaparsecs, the universe shows an in- 
tricate web of filaments and walls composed of densely clus- 
tered galaxies which are surrounded by comparatively empty 
voids. These voids appear to be the largest coherent structures 
in the universe. Their size, though vast enough in itself, is 
still but a small fraction of the observable universe today. It 
is therefore often thought that inhomogeneous structures such 
as these voids should have little influence on the overall ex- 
pansion of the universe. 

One reason to question this assumption, or at least to wish 
to test it through specific instances, is that the equations of 
general relativity governing the mutual influence of matter on 
the shape of space-time are highly complicated and inherently 
nonlinear. Although the universe has a few naturally avail- 
able small parameters — the initial smallness, for example, of 
the fluctuations in the density of the matter, or the ratio be- 
tween the typical size of a void and the size of the observable 
universe — it is also a dynamic and an evolving system. The 
uniform part, thought to determine the overall expansion of 
the universe, and the spatially varying part, which describes 
the formation of structures such as the voids and the clusters 
of galaxies, grow at very different rates. Even if the average 
effect of the inhomogeneities in the universe is extremely tiny 
at an early time, this smallness can be undone once the dif- 
ference between how these two parts evolve is given enough 
time to show itself. 

This article will demonstrate how the growth of the inho- 
mogeneities in a universe containing only matter can affect 
its average expansion rate. The example that will be used to 
illustrate this behavior will not be anything quite as compli- 
cated as the irregular system of walls and voids that occurs in 
the actual universe, but will instead consist of a universe com- 
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posed of matter that has a single periodic fluctuation with a 
fixed comoving wavelength. Also, the fluctuation will be cho- 
sen to depend on only one of the spatial directions. With this 
simplicity, the growth of the fluctuation in this example can be 
solved analytically, at least until the time when the distribution 
of the matter becomes markedly inhomogeneous. Beyond this 
point, the analytic solution becomes a poor approximation to 
the true behavior. But despite its simplicity, this example is 
meant to capture several properties of the actual universe. It 
begins in a nearly uniform state and evolves into a regular 
pattern of thin walls and wide voids as the matter is drawn 
toward the initially slightly denser regions and away from the 
less dense places between. 

In this simplified version of the universe, it is possible to 
solve for the detailed time-dependence of the inhomogeneous 
part well beyond a linear approximation, and its periodicity 
means that it is straightforward to define its spatially aver- 
aged behavior. An important result of this analysis is that it 
provides a way to characterize when the expansion of the uni- 
verse begins to depart strongly from the behavior expected 
of one composed of uniformly distributed matter. The con- 
dition that separates these phases can be stated in terms of 
two dimensionless parameters. If the initial amplitude of the 
spatial variations in the density of the matter, as a fraction 
of its average density, is called 8, and if the size of the ob- 
servable universe measured in units of the wavelength of the 
fluctuation is called n, then the advent of the stage when the 
inhomogeneities strongly influence the overall expansion rate 
occurs when n 2 8 ~ 1. The average expansion rate during the 
latter stage is faster than the earlier one, so as this transition 
occurs between them, the average expansion rate of the uni- 
verse begins to speed up. Using the measured size of the ini- 
tial fluctuations in the density of the matter, if this transition 
were to have occurred approximately four billion years ago, 
the corresponding wavelength of the fluctuations today would 
be of the same size as the voids in the universe. Whether this 
acceleration in the average expansion rate actually translates 
into a luminosity distance function that resembles that of the 
observed universe will be treated in later work, once the exam- 
ple here has been generalized slightly to be made to resemble 
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more closely the observed universe. 



II. A SINGLE FLUCTUATION 

From what has been observed of the universe so far, it 
appears to have begun in a largely homogeneous state that 
evolved into a truly inhomogeneous one only long after all 
the other known substances within it had been so diluted, in 
comparison with the density of the matter, that they had only a 
negligible role in its later expansion. It is therefore a good ap- 
proximation to neglect from the start the radiation — and any 
other possible ingredients — and to treat the universe as though 
it were composed entirely of nonrelativistic matter, which is in 
the form of a pressureless idealized gas. This gas will not be 
perfectly uniform throughout space, but will have some spa- 
tially dependent fluctuations even at quite early times. To sim- 
plify the analysis a little, this article will consider a case where 
the inhomogeneities exist in only one of the spatial directions. 
Choosing a set of coordinates that is consistent with these con- 
ditions, this space-time can be described by the metric, 

ds 2 — g^ v (t,x) dx^ dx v 

= c 2 {t 1 x)dt 2 -b 2 {t,x)dx 2 -a 2 (t,x) [dy 2 + dz 2 ]. (2.1) 

The important property of the matter for its effect on the 
space-time is its density, p(t,x). The energy and the momen- 
tum associated with this pressureless gas is 
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c 2 (t,x)p(t,x)dt 2 , 
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and their conservation corresponds to the condition V^T^ = 
0. This equation imposes two conditions on the time-evolution 
of the density and on the spatial dependence of the time com- 
ponent of the metric, 
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and 
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The integrated form of the first equation becomes the state- 
ment that the only time-evolution in the matter's density is 
that produced by its dilution through the expansion of space, 
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The universe is never completely empty, p ^ 0, so the sec- 
ond of the conservation conditions indicates that the function 
c(t ,x) — > c(t) is independent of where it is being evaluated. It 
can then be dispensed with altogether by redefining the time- 
coordinate to absorb this function, c(t)dt — > dt, leaving just 

ds 2 =dt 2 -b 2 (t,x)dx 2 -a 2 {t,x) [dy 2 + dz 2 ]. (2.5) 

How the geometry of the space-time responds to the pres- 
ence of the matter is described by Einstein's equation, 



The components of this equation translate into three separate 
conditions on the evolution and the spatial dependence of the 
metric functions, a(t ,x) and b(t,x), 
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Now to go further, it is necessary to choose a particular form 
for the comoving density, po(x), and the early behavior of the 
metric. 

The approach to be followed here is complementary to the 
standard treatment of the growth of the universe from a highly 
uniform state into a very inhomogeneous one. Initially, the 
spatially dependent fluctuation in the metric and in the mat- 
ter's density are both extremely small in comparison with the 
uniform parts of each. In this state, the universe can be ap- 
proximately described as a small inhomogeneous piece that 
evolves in a homogeneous background, which is often written 
in a standard form, 



ds 2 = dt 2 -r 2 {t) [dx 2 + dy 2 + dz 2 



(2.8) 



It is then assumed that the growth of the inhomogeneous 
component — as the matter begins to accumulate into regions 
of higher density leaving emptier voids among them — while 
important at scales comparable to the size of the fluctuations, 
has a negligible effect on the expansion rate of the universe on 
much larger scales. 

Yet the actual behavior might be quite different. The rela- 
tivistic dynamics governing both the overall expansion and the 
collapse of the matter are highly nonlinear and include parts 
evolving at different rates. So while the universe has a scale 
which is small, the initial smallness of the fluctuations, there 
are other time-dependent factors, such as the relative growth 
rates of the homogeneous and inhomogeneous parts of the 
metric. After a sufficient time, the latter can overwhelm the 
former. 

So rather than following the growth of a full spectrum of 
fluctuations in a fixed background, this article will instead fol- 
low the evolution of a fluctuation with a single wavelength 
without separating it from the 'background.' This example 
will show directly how and when the growth of the inhomo- 
geneities alters the expansion of the universe at comparatively 
large scales. 

With this picture in mind, consider the evolution of what 
would initially appear to be a sinusoidal fluctuation in the den- 
sity of the matter, 



PoO) =po[l + ecos(2jrfar)], 



(2.9) 



(2.6) 



which would be a single Fourier mode of a more general spec- 
trum of fluctuations. The average comoving density of the 
matter is po- This fluctuation is meant to evolve into some 
of the features of the large-scale structure of the universe to- 
day. However, the universe should have been largely uniform 
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over all scales when, for example, the physical wavelength of 
the fluctuation was just equal to the size of the Hubble hori- 
zon. So, in the beginning, the amplitude £ provides a good 
parameter for describing the metric perturbatively; to a first 
approximation the components of the metric do not depend 
on the spatial coordinates at all, and are moreover equal to 
each other, 



a(t,x) = d t 2/3 + 0(e) 
b(t,x) = d t 2/3 + 0(e), 



(2.10) 



where al — 6nGpo- By the Hubble horizon, what is meant is 



the horizon associated with this homogeneous part, defined by 
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Although this scale will also be later used as a yardstick, 
once the inhomogeneities have grown sufficiently large it be- 
comes increasingly less meaningful as a real, physically rele- 
vant scale. 

Having this small parameter available permits the metric to 
be calculated analytically. Up to corrections proportional to 
e 5 , the solution to Einstein's equation is 
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where the terms with the same time dependence have been 
grouped together. The inhomogeneous terms are obviously 



growing at a faster rate, t 



4/3 f 2 



than the homogeneous 



one, f 2 ' 3 . So while they start quite small, being suppressed by 
various powers of £, with enough time they eventually become 
the more important part for the expansion of the universe. 

To gain a clearer insight into the precise time at which the 
nonlinear, inhomogeneous regime begins, the time coordinate 
can be replaced with a more natural unit. The comoving wave- 
length of the fluctuation provides one natural scale in this uni- 
verse; a second one is given by the size of the comoving (ho- 
mogeneous) Hubble horizon, dn/a \ e =o- Measuring the latter 
in units of the former produces a dimensionless ratio, 
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At times, it is a little more convenient to absorb certain recur- 
ring factors, so define N(t) to be a bit more than four times 
this ratio, 



N(t) 
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Counting the number of wavelengths per horizon thus, the 
leading behavior for each term of a distinct time-dependence 



in the metric becomes 

a(t ,x) = aof 2/3 { 1 + sin 2 (27Tfcc) 

-2ro^ 4 e 4 sin 4 (2^x) + ---| 
b(t,x) =fl f 2//3 jl - y^V 2 ecos(27z;fc<:) 

+ ^ 7Y 4 e 3 cos (2 nkx) sin 2 (2nkx) + • • • | . (2.15) 

The departure from approximate homogeneity occurs first in 
the x-direction, as should have been anticipated, when the 
number of wavelengths per horizon reaches N <~ e -1 / 2 . 

in. A POSSIBLE MISATTRIBUTION 

The question now arises — what mistake is made in assum- 
ing that the growth of the inhomogeneities has no effect on 
the large-scale expansion of the universe? While the inho- 
mogeneous terms evidently grow much faster than the homo- 
geneous part, it might be thought that on average their effect 
vanishes — if the extra density in one place tends to slow the 
expansion a little, the surrounding voids will hasten it. At 
linear order in the amplitude e, since all of the spatial depen- 
dence is a cosine, the spatial average of this function does 
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vanish; however, this property is not true for the higher order 
terms. Even after averaging over the spatial coordinates, the 
effects of the inhomogeneities will remain. If one then insists 
on seeing the large-scale expansion as that produced by uni- 
formly distributed matter, it will appear that the expansion has 
unaccountably begun to speed up. 

The periodicity of this fluctuation means that it is quite sim- 
ple to define spatially averaged quantities. For a general func- 
tion f(t,x), define /(f) to be the result of averaging it over a 
wavelength, 



Since the metric has a spatial dependence on only one of 
the coordinates, it will remain anisotropic even after aver- 
aging over the x-dependence. Defining therefore two scal- 
ing functions a(t) and b(t) for the averaged metric, g^ v {t) = 
(8nv(t,x)), 



ds 2 = dt 2 - b 2 (t)dx 2 -d 2 {t) [dy 2 + dz 2 



(3.2) 
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Among the terms with the same time-dependence, which is 
implicit in the N(t) dependence, the leading terms at a partic- 
ular order in N and e are 
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Although a(f ) has no true analogue in the actual universe, 
since there is no genuine translation invariance in any direc- 
tion, even its behavior is seen to stray from a f 2 / 3 expansion 
when N grows to be the same size as £ _1 . The importance of 
the inhomogeneities appears much sooner in the x direction, 
when N ~ e -1 / 2 . When comparing to physical scales in the 
universe, it is important to remember that 7Y(f) is not quite 
equal to the number of wavelengths per horizon, which was 
instead called «(f), nor is e equal to the initial amplitude in 
the physical density of the matter. If at the time fo, the wave- 
length of the fluctuation is exactly equal to the Hubble horizon 
(n(fo) = 1), then the amplitude of the initial fluctuation in the 
physical density, 
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which has been called 8 to distinguish it from the amplitude 
in the comoving density, is a bit larger than e, 
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Replacing N and e with n and 8 does not much alter the form 
of the scaling factor, 

b(t) = a f 2 / 3 {l + \n A 8 2 - ±n H 8 4 + •••}. (3.7) 

Once n 2 8 has grown larger than about 2 /3, the expansion rate 
begins to increase. 



Remember that the horizon size tends to grow at a faster 
rate than the physical wavelength of the fluctuation, at least 
during the early stage when the expansion rate is governed by 
the homogeneous term. Beyond this stage, while n(t ) can still 
be formally defined, it corresponds less and less to the mean- 
ing previously assigned to it. But since this departure sets in 
fairly gradually, and since from a cosmological perspective it 
began comparatively recently, n should still very roughly cor- 
respond to the number of wavelengths per horizon. 

The threshold between the era governed by the homoge- 
neous term and that when inhomogeneities become more im- 
portant for the overall expansion occurs when n 2 8 <~ 1. Just 
prior to this threshold, the expansion rate begins to accelerate 
as it passes from a t 2 / 3 dependence to a f 2//3 n 4 <~ f 2 depen- 
dence. Beyond this point, the solution here is not sufficient 
for determining the extent of this acceleration and whether 
and how long it lasts, since a perturbative description of the 
metric based solely on the smallness of 8 or e fails. 

To see, from a slightly different perspective, how the depar- 
ture from homogeneity leads to an acceleration in the expan- 
sion rate of the universe, consider the first of the components 
of Einstein's equation, which in a homogeneous universe cor- 
responds to Friedmann's equation, 



Goo = ^oo - 3^oo^ = 87tG Too. 



(3.8) 



The spatial average of the left side leads to a more compli- 
cated apparent expansion than the f~ 2 scaling that could be 
explained by a uniformly distributed gas of matter alone, 
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Most of these terms are only tiny corrections, or become im- 
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portant only far into the regime when the inhomogeneities 
have grown large. Immediately before the transition into this 
regime, the leading part of the gravitational part of the aver- 
aged Friedmann's equation is 

Goo(0 = lji{ 1 + \ n " 82 + I"* 5 * + •••}• < 3 - 10 ) 

If one insists on accounting for this behavior only through uni- 
formly distributed materials, ordinary matter alone would ap- 
pear not to be sufficient. By itself, a uniform gas of matter 
would explain the first term, but other materials would need 
to be invoked to account for the other terms, both of which 
are diluted more slowly than the first. The second term scales 
as n 4 /f 2 ~ f~ 2 / 3 , while if some strange substance were in- 
troduced to explain the third, its density would appear to be 
growing with the expansion, as n & /t 2 ~ f 2 / 3 . However, this 
would be a little misleading; once the universe has passed the 
threshold were n 2 8 ~ 1, the actual behavior of this artificial 
substance would be determined by adding together all of the 
terms scaling as arbitrary powers of n 2 8, at least until the next 
threshold when the terms scaling as powers of n 3 5 2 become 
important too. 

In this example no such material is needed. It would only be 
invoked by insisting that the large-scale expansion is governed 
by a uniformly distributed substance and thereby misattribut- 
ing the origin of the change in the expansion rate. The source 
of the acceleration and the advent of the new expansion rate 
are due here to the inhomogeneities in the distribution of the 
matter alone. 

The fact that the average of the scale factor accelerates does 
not by itself immediately mean that the luminosity distance 
function, which is measured through the observations of dis- 
tant supernovae, behaves as it does in the actual universe. The 
calculation of this function will be left to later work, since 
it is important first to understand the behavior well into the 
inhomogeneous era, which probably needs to be studied nu- 
merically, to see the extent of the acceleration, both in terms 
of its overall magnitude and its duration. 



IV. PHYSICAL SCALES 

This example of an inhomogeneous universe was con- 
structed to be a simple as possible, while still keeping the es- 
sential features of the actual universe that were to be analyzed. 
The results of the calculation showed that once the inhomo- 
geneities have matured sufficiently, they can significantly af- 
fect the large-scale expansion, which begins to accelerate as 
it passes into this later stage. This effect is found by calcu- 
lating the higher order, nonlinear behavior of the solution to 
Einstein's equations; but because of the inherent smallness of 
the initial amplitude of the fluctuations, a fair amount of time 
must first pass before this influence appears, during which the 
horizon grows much larger than the wavelength of the fluc- 
tuation. So the idea that because the resulting structures are 
small compared with the size of the observable universe they 
are unimportant is a mistaken one — at least in this example — 



or rather it is true only while nonlinear parts of the solution 
are still small. 

The expansion of the actual universe seems to have de- 
parted from the behavior expected of a space-time uniformly 
filled with matter about four billion years ago or so. If this 
departure was produced by some of the regularly occurring 
larger structures in the matter, then their size at that time 
should be simply related to the amplitude of the initial fluc- 
tuation, 8. From what can be inferred from the observa- 
tions of the Wilkinson Microwave Anisotropy Probe [|l|], the 
size of the product of the change in the spatial curvature 
from its average value, evaluated at two different places, is 
(2.430 ± 0.091) x 10~ 9 , with only a slight dependence on the 
distance between the two points. If the initial amplitude of 
the matter's density is taken to be comparable to the square 
root of this value, so that 8 ~ 5 x 10~ 5 , then the stage of 
the expansion that is strongly affected by the inhomogeneities 
ought to have begun when the structures responsible for it had 
a typical size so that n ~ 140. To evolve their size from its 
value four billion years ago to today requires knowing the 
time-dependence of the metric in the era governed by the in- 
homogeneities, which is not possible for the approximation 
used in this calculation. However, assuming that the expan- 
sion rate, which has only changed fairly recently according to 
a cosmological reckoning, does not depart excessively from 
the n ~ t '/ 3 dependence of the homogeneous era before it, the 
size of the structures responsible for the acceleration would 
now be '/i60 the size of the horizon. The appropriate hori- 
zon is not the size of the observable universe today, which 
is estimated to be about 14.5 Gpc, but is rather the Hubble 
horizon. This horizon is inversely proportional to the current 
value of the Hubble parameter, Ho, used to compute the rate of 
the recessional velocity of distant objects. It is also measured 
under the assumption of a homogeneous cosmology on aver- 
age. The size of this Hubble horizon, c/Hq, is approximately 
4 Gpc, which in turn implies that the size of the structures that 
are chiefly responsible for the change in the recent expansion 
history of the universe ought to be about 30 Mpc today, which 
is comparable to the size of the smaller voids. 

In the actual universe, the voids are fully three dimensional 
and not periodic in one dimension only, nor are they arranged 
in a simple periodic array, nor is there a unique size of the 
voids. Their size today typically falls within the range of 
about 40-90 Mpc [2]. When the Hubble horizon is measured 
in units of this range of void-widths, and then followed back 
to four billion years ago, still very crudely applying the t 1//3 
dependence, the range for the values of n falls between 40 and 
90. This range corresponds to a range of values for n 2 8, 

0.1 <n 2 8 <0.4, (4.1) 

which is consistent with the expectation that the acceleration 
should have begun once n 2 8 ~ 2 /3. 

Whether or not this mechanism is in fact responsible for 
the relation that has been observed between the luminosities 
and the redshifts of distant supernovae requires improving the 
example presented here so that it includes a few more of the 
properties of the actual universe. Enlarging the spatial de- 
pendence to include all of the coordinates and allowing for a 
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less regular array of voids probably will not alter much the 
conclusions from this simpler periodic example. It is proba- 
bly more important to include a full spectrum of fluctuations. 
Shorter wavelengths have been growing for a longer time than 
the scales associated with the voids, so it might be expected 
that the effects of the inhomogeneities at smaller scales ought 
to have appeared earlier. In its full, three-dimensional mo- 
mentum representation, the power spectrum for the primordial 
fluctuations inferred from observations decays approximately 
as k~ 3 . A smaller amplitude for shorter wavelengths would 
provide a natural suppression in 5 that would compensate for 
the longer time, n, that shorter wavelengths would have been 
smaller than the Hubble horizon. This effect would addition- 
ally reduce the role of mixings among different wavelengths, 
which would also need to be evaluated when the metric is 
treated well beyond a linear approximation. Once these prop- 
erties have each been addressed, it should then be possible to 



evaluate how the true inhomogeneities affect what is actually 
observed. 

The encouraging conclusion of this work is that it shows 
that it might not be necessary to abandon altogether the idea 
that the universe is composed primarily of matter nor that 
its dynamics are governed by ordinary general relativity, in 
which case the discrepancies between some of the earlier the- 
oretical expectations and observations would rather be due to 
having insisted on assumptions that do not hold in the actual 
universe. 
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